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ABSTRACT 

One can define a random walk on a hypercubic lattice in a space of integer dimension 
D. For such a process formulas can be derived that express the probability of certain 
events, such as the chance of returning to the origin after a given number of time steps. 
These formulas are physically meaningful for integer values of D. However, these formulas 
are unacceptable as probabilities when continued to noninteger D because they give values 
that can be greater than 1 or less than 0. In this paper we propose a random walk which 
gives acceptable probabilities for all real values of D. This D-dimensional random walk 
is defined on a rotationally-symmetric geometry consisting of concentric spheres. We give 
the exact result for the probability of returning to the origin for all values of D in terms 
of the Riemann zeta function. This result has a number-theoretic interpretation. 
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I. INTRODUCTION 



Recent studies of quantum field theory in D-dimensional Euclidean space have fo- 
cussed on the possibility of expanding the Green's functions of the theory as series in 
powers of Z).^'^'^ A technique used in this work is to introduce a Z)-dimensional hyper- 
cubic lattice, and, using strong-coupling lattice methods, to obtain expressions for the 
Green's functions for integer values of D. Order by order in these strong-coupling ex- 
pansions the coefficients are polynomials in D and therefore appear to have simple and 
natural continuations to noninteger values of D. Of course, such a procedure is not unique; 
there may be many different analytic continuations off the integer-dimensional hypercubic 
lattices. The procedure was justified on the basis of comparisons with known results in 
solvable quantum field theoretic models. 

Such an analytic continuation does not always produce acceptable physical results. For 
example, if we try to define a random walk in arbitrary-dimensional space by analytically 
continuing a random walk off an integer-dimensional hypercubic lattice, we can obtain 
unphysical values for probabilities for certain physical processes. To illustrate very simply 
how such problems can arise, consider the probability P2 that a random walker on a lattice 
in L>-dimensional space will return to the starting point after two steps. To compute this 
probability we consider D to be a positive integer and the lattice to be hypercubic. For 
such a lattice, 

P2 = — ■ (1-1) 

^ 2D ^ ' 

The obvious analytic continuation of this result to noninteger D, namely the function 
1/(2D), cannot be the correct probability when D < 1/2. 

The purpose of this paper is to devise a D-dimensional geometry in which a random 
walk is well defined when D is noninteger. In this model all probabilities have sensible 
values; they lie between and 1.^ 

This paper is organized as follows: In Sec. II we review the conventional treatment of 
random walks on integer-dimensional hypercubic lattices. We show that if we attempt to 
continue the results to noninteger dimension we obtain unacceptable values of probabilities 
for physical events. In Sec. Ill we invent a spherically symmetric geometry on which a 
random walk is well defined for all real values of D. This geometry eliminates the need 
for analytic continuation. We show how to express the probabilities in this model as 
solutions to a difference equation. We solve this equation and show that for all real D, 
probabilities have physically acceptable values. One quantity of particular interest is Hi), 
the probability that a random walker will eventually return to the origin. We determine 
Ho in closed form for all values of D. For one interesting special case, in which the random 
walk occurs on a geometry of concentric spheres whose radii are consecutive integers, 

"^"^"OD^ (^>2) , (1.2) 

where C is the Riemann Zeta function. This result has an interesting number-theoretic 
interpretation. 
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II. RANDOM WALKS ON HYPERCUBIC LATTICES 



In this section we review the standard formalism^ for describing random walks on a 
D-dimensional hypercubic lattice, where D is an integer. Let Cn,t represent the probability 
of the random walker being at the position n at time t. Note that n is a D-component 
vector having integer entries and that t is a nonnegative integer. The random walk is 
assumed to begin at the origin at t = so that 

Cn,o = (^n,o • (2.1) 

The random walk consists of unit-length steps connecting adjacent vertices on the lattice. 
Thus, at each lattice point the walker may go in any one of 2D directions with equal 
probability. 

To describe a random walk we construct a recursion relation satisfied by the proba- 
bilities Cn,t: 

1 ^ 

<^n,t = ^ ^ (Cn-ej,t-l + C'n+ej^t-l) , (2-2) 

where e-,- is the Z)-component unit vector for which all entries are zero except for the jth. 
component which is 1. This equation states that the chance of being at the position n at 
time t is the sum of the probabilities of having been at 2D equally- likely neighboring sites. 
To solve this recursion relation we introduce the discrete Fourier transform; to wit, 

I(Cf'^*(^)^"'''" ' 

where k is now a continuous D-dimensional vector. In terms of A't(k) the recursion relation 
reads 

1 ^ 

A't(k) = A'i_i(k)-5^cosA;,- . (2.4) 



Cn,t — 



D 



The general solution to this equation is 



A',(k)=A(k) l^if^cos^j . (2.5) 

To determine the constant ^(k) we use the initial condition in (2.1) which specifies that 
^(k) = 1. We now obtain Cn,t by taking the inverse Fourier transform 

t 

„— ik-n 



= Sin ^"4:1 



x=0 
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where I,^{x) is the modified Bessel function. 

We are particularly interested in the probability of returning to the origin n = 0: 



2t 

(2.7) 



where 



|2t . ■ ' ^ 

X 



1=0 ^ ' 

(Note that Co,2t+i = 0.) For D = 1 for instance, we obtain the familiar result 

Co,2t = , (2.8) 

for the probability of a one- dimensional random walker to return to the origin after t time 
steps. 

For arbitrary integer D we have 

U,2 - ^ , 

Co,4 = ^(6^-3) , (2-9) 
^0,6 = ^ (6i;^ - + 4) . 

Observe that for some noninteger values of D these expressions are not acceptable probabil- 
ities. Specifically, Co,2t is oscillatory until D is sufficiently large, at which point it is positive 
and monotone decreasing. On Fig. 1 we have plotted Co, 2, Co, 4, and Co,6- These values 
become acceptable as probabilities only when D > 1/2, D > 1/2, and D > 0.7224..., 
respectively. In general, as t increases, the lowest value of D at which Co,2t becomes an 
acceptable probability also increases. On Fig. 2 we plot scaled versions of Co,i8 and Co, 20 
as functions of D. Note that Cq.is becomes an acceptable probability only when D is 
roughly 2 or larger, and Co,2o only when D is roughly 2.6 or larger. A more precise picture 
of the behavior of Co,2t is given in Table 1. In this table we give the critical values of the 
dimension ^critical for which the function Co,2t can first be considered a probability. That 
is, for aU D > L>criticai, < Co,2t < 1- 

Having constructed the functions Co.2f we can now construct expressions P2t which 
represent the probability that the random walker returns to the origin at time 2t for the 
first time. The first few such expressions are 

P2 = Cq 2 , 
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P4 = Co 4 — Cq 2 , 
Pq = Cqq — 2(70,2^^0,4 + 1 
-Ps = C'0,8 — 2Co, 2^^0,6 — 4 + 3Cq 2(70,4 — ^^0,2 • 

In general the relation between P2t and (7o,2t is given by the elegant formula^ 

00 ^ 

Y.P^t^' = ^- ^oo ^ . . (2.10) 

Also, if we evaluate the above generating functions at a; = 1 we obtain the probability IId 
that the random walker ever returns to the origin: 

oo ^ 

nz, = ^P2t = 1- ^00 ^ • (2.11) 

The quantity in (2.11) represents a physically interesting probability, namely the 
chance that a random walker will eventually intersect her starting point. However, as we 
can see, its mathematical form is not very satisfying because it is presented in terms of a 
sum over probabilities (7o,2fc) each of which is an unacceptable probability for some range 
of noninteger D. Nevertheless, we can perform a formal Borel-like summation of the series 
in (2.11) by substituting the integral representation for (7o,2t in (2.7) and interchanging the 
orders of summation and integration. We obtain immediately a closed-form quadrature 
expression for II d: 

We emphasize that the expression in (2.12) for Hd derives from an infinite sum over 
probabilities that are unacceptable. However, when recast in the form in (2.12) we find 
the surprising result that H^j is itself an acceptable probability for all real positive D if 
we interpret the expression as follows. For < D < 2 it is clear that the integral in (2.12) 
diverges. Thus, we conclude that 

Ud = 1 (0 < L> < 2) . (2.13) 

This conclusion verifies the well-known result, sometimes referred to as Polya's theorem, 
that a random walker in one or two dimensions will ultimately return to the origin with a 
probability of 1. 

For D > 2 the integral in (2.12) exists and is a smooth function of D. We can thus 
regard (2.12) as a deGnition of H/^ for noninteger D. We find that for L> > 2, Hd is a 
smooth, monotone-decreasing function of D and is thus an acceptable probability function. 
For example, for the special case D = 3 we obtain^ 

H3 = 1 — = 0.3405373296 ... . (2.14) 

v^r(l/24)r(5/24)r(7/24)r(ll/24) 
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Unfortunately there is no way to evaluate the integral in (2.12) for arbitrary D or 
even arbitrary integer D. Thus, we are forced to study this integral representation asymp- 
totically. For large D, Hd has the asymptotic form 



1 1 7 35 215 1501 5677 44489 



{D ^oo) . (2.15) 



2D 2D'^ 8D^ leD"^ 32D^ 6AD^ GAD^ 128D» 



This expansion is not accurate when D is as small as 3. However, it rapidly becomes 
accurate as D increases. The structure of (2.15) is not as simple as it first appears; 
eventually the signs of the terms become negative, remain negative for many orders, become 
positive for many orders, and appear to continue oscillating in this block pattern. 

Note that from the above description, the function II^, is also not smooth at D = 2; 
it has a discontinuous derivative there. When D < 2, Hd = 1. On the other hand, when 
D>2 



Thus, the function in (2.12) approaches its value at -D = 2 from above with the negative 
slope — 7r/2. 

If we analytically continue the function in (2.12) from values of D above 2 to values 
below 2 we obtain an unacceptable probability; wc find that Ho rises as D decreases and 
eventually becomes infinite at D = 1 (see Appendix). All of the above properties of the 
function Hd are summarized in Fig. 3. 

Our objective in this paper is to construct a lattice geometry for which the function 
Hd is an infinite sum of terms, all of which are acceptable probabilities for all real D and 
not just integer D. Our construction is given in the next section. 




-{D-2) + ... {D^2+) 
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III. RANDOM WALKS IN D DIMENSIONS 



To define a random walk in D dimensions, where D is not necessarily an integer, we 
return to the discussion of a random walk on a hypercubic lattice. Normally, we regard a 
random walker as going from one lattice point on this lattice to an adjacent point. However, 
for a completely equivalent way of viewing the walk, consider the dual lattice in which each 
lattice point is replaced by a D-dimensional hypercube of side 1. The walker begins in the 
center cube; we define region 1 as the interior of this cube. There are 2D cubes contiguous 
to region 1 into which the walker may go on the first step. The region inside these 2D 
cubes is labeled 2. On the second step the walker may return to region 1 or she may enter 
region 3, the region composed of all cubes contiguous to region 2 except for region 1. 
Region 3 comprises 2D^ cubes. On the third step the walker may return to region 2 or 
she may enter region 4, the region composed of all cubes adjacent to the cubes in region 3 
except for those cubes in region 2. There are 2D{2D'^ + l)/3 cubes in region 4. In Fig. 4 
we depict the first five of these regions for the case D = 2. Note that if the walker is in 
region n, she must enter either region n + 1 or region n — 1 on her next step. At each time 
step the walker measures the area of the surfaces bounding the hypercube she occupies 
and travels outward or inward with a probability in exact proportion to the outward or 
inward surface area she must traverse. 

It is easy to visualize this model of a random walk in any geometry. The simplest 
such geometry is, of course, a set of regions bounded by concentric nested spheres. As we 
will see, such a geometry docs not single out the integer dimensions as being special in any 
way; a random walk can be defined in any real positive dimension D.^ 

Specifically, we consider an infinite set of concentric nested spheres of radius Rn 
(n =1, 2, 3, . . .). Then region n is the volume between the n — 1st and the nth spheres 
(see Fig. 5). Thus, if the walker is in region n the probability of walking outward, Pout('T')) 
and the probability of walking inward, Pin{n), are given by 

Pout(n) - ^^(^-^ 



SoiRn) + SD{Rn-l) ' ^2 -j^^ 

Pi„(n) = SpjRn-l) 

SD{Rn) + SniRn-l) 

where Sd{R) is the surface area of a D-dimensional sphere of radius R: 

'^^^^^ = 1X^72)^''"' • 

Since the surface area <S'o(-R) is proportional to R^~^, the inward- walk and outward- walk 
probabilities are 

rD-1 

Pout{n) = p-D-l , r>D-l ' 

r,D-l (3-2) 

p. (n) - 



7 



These formulas are valid for all D and all n > 2. For n = 1, Pout(l) = 1; and -Pin(l) = 
because there is no region inward from n = 1 and therefore the walker must go outward. 
Note that Pout and Pin are nonnegative numbers less than or equal to 1 with Pin(n) + 
Pout{n) — 1; thus, they are acceptable probabilities for all real D7 

In this model, the probability Cn^t thai the random walker arrives at regionn at time t 
satisfies a recursion relation analogous to that in (2.2), 

Cn,t ^ Pin{n + l)Cn+l,t-l + Poutin - l)Cn-l,t-l (n > 2) , ^^^^ 
Cl,t = -Pin(2)C2,t , 

where Pin and Pout are given by (3.2). This equation states that for a random walker to 
be in region n at time t, she must have come from region n + 1 or region n — 1 at time 
t — 1. The initial condition satisfied by Cn,t is 

Cn,0 = Sn,l , (3.4) 

which states that at time t = the random walker is located in the central sphere 
{region 1). 

We can now use (3.3) and (3.4) to calculate Ci^2, the probability of being in the central 
sphere {region 1) at the second time step: 

p-D-l 

^1,2 — 



p-D-l I D-D-1 



Since P2 > -Ri > we conclude that < Ci^2 < 1 and thus that Ci 2 [unlike Co, 2 in 
(2.9)] is an acceptable probability for all real values of D. For the special simple case of 
equally-spaced concentric spheres for which P„ = n we have plotted Ci,2, ^1,4, Ci,6 in 
Fig. 6. Note that these probabilities do not oscillate like those in Figs. 1 and 2 and they 
are acceptable probabilities for all real D. 

We present two methods for obtaining from (3.3) the value of IId, the probability 
that a random walker who starts at the central sphere eventually returns to the central 
sphere {region 1). The first technique involves converting the partial difference equation 
in (3.3) to an ordinary difference equation by summing over t. We begin by making the 
substitution 

Cn,t = {R^-' + R^ll)Er,,t (n>2) , 
Ci,t = Ri ^Ei^t • 



(3.5) 



Then the equation satisfied by t is 



(P„ + Rn-l )En,t — Rn -E'n+l,t-l + Rn-1 En-l,t-l («■ > 2) 

Ei,t = E2,t-i , 



(3.6) 



where we have substituted the values of Pin and Pout in (3.2). In addition, the first line in 
(3.6) is also valid for n = 1 and D > 1 if we define Rq = 0. 
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We can eliminate the t index from (3.6) by summing over t. We define 

oo 

Fn = J2En,t (n>l) (3.7) 



t=0 



and sum (3.6) from t = 1 to t = oo. [Of course, we must assume here that the result is 
finite. If the sum over Cn,t is infinite then from (2.11) we have immediately that Hd = 1-] 
The result is 



(3.8) 



(R^ ^ + R^_l)Fn — Rn + R^_l F^-i (n > 2) , 

i?f-i (Fi - F2) = 1 , 

where we have made use of the initial condition in (3.4). 

The first line in (3.8) can be rewritten in exact-diff'erence form: 

[R^-'F^ - = - R^llFr,] (n > 2) , (3.9) 

which we then sum from n = 2ton = A'" — Ito give 

R^ZIFn-i-R?-'Fi^R^ZIFn-R?-'F2 (iV>3). (3.10) 

We eliminate F2 from (3.10) using the second line in (3.8) and obtain 

[Fn-Fn-i] = -R^^Fi~^^ (iV>2) (3.11) 
and sum it from N = 2 to N = K. The result is 

Fk = F,-J^-^ • (3-12) 

N=l 

If Fk ^ as -ftT — > 00 (we will establish the validity of this limit later) then we have 
the simple result that 

00 cx) ^ J-, X D—1 

EC-u-E^ (3.13) 

t=0 N=l ^ ^ 

from which we can express the probability that a random walker returns to the origin as 
a generalized zeta function 

^D = l ^--d3T • (3.14) 



More generally, for arbitrary functions Pin{n) and Pout(^) with Pin(n) + -Pout('^) = 1; 
n > 2, and Pin(l) = and -Pout(l) = 1, we obtain from (3.3) 

00 00 Af j-j / N 

t=0 iV=2 n=2 ^outl.'^t; 



and 

IId = 1 ^7-^ ■ (3-16) 

-L + Z^JV=2 lin=2 Pout(n) 

The results in (3.15) and (3.16) reduce to (3.13) and (3.14), respectively, for the choice of 
Pi„(n) and Pout{n) in (3.2). 



Alternative Solution in Terms of Continued Fractions 

If we simply iterate (3.3) we obtain a sequence of increasingly complicated formulas 
for Ci^2t- 

Ci,o = 1 , 

' f3 17) 

Ci,4 = Qi(gi + g2) , ^ ■ ^ 

Ci,6 = Ql[(gi+Q2)'+Q2Q3] , 

and so on, where 

Qn = Pinin + l)Poutin) . (3.18) 

However, the sequence of formulas in (3.17) is a clearly recognizable pattern:^ Ci^2t are the 
coefficients in the Taylor series for the function whose continued-fraction representation 
has the coefficients Qn- 



oo 
t=0 



Thus, 



and 



Ci,2tx' = 1/(1 - Qia;/(1 - Q2x/{1 - Q^x/{\ - Q^x/ . . .)))) . (3.19) 



= 1/(1 - Qi/{1 - Q2/{1 - Qs/il -Qa/.. .)))) (3.20) 

t=Q 



liD = Qi/(1 - Q2/(i - Qs/ii - Q4/(i - Qs/ • • ■)))) • (3.21) 



From (3.20) it is easy to recover (3.15) for arbitrary Pin and Pout, or (3.13) using 

Q^ = Pf -V(i?f + R?-') and Qn = Rl^'-ViiR^ll + R^-'){R^-^ + <+7)], ^ > 2, 
by the following procedure. For instance, if we truncate (3.20) before the first minus 
sign we obtain the first term in the sum on the right side of (3.13); if we truncate before 
the second minus sign we obtain the first two terms in the sum; and so on. Thus, as a 
byproduct of this analysis we have established the continued-fraction identity 

y^^ = i/ f 1 - ^1 /(i ^1 /(i ^1 /... 

^^Sn '\ s^ + sJk {s^ + S2){S2 + s^)' \ {S2 + s^){s^ + s^y 

= 1/ (1 - S^/ {Si + 52 - si/ {S2 + ^3 - si/ (^3 + -54 - si/ ...)))) , 

(3.22) 
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which contains as a special case a continued-fraction representation of the zeta function: 
C(r) = 1/(1 -1/(1 + 2'" -2^7(2^ + 3^ -3^7 (3'" + 4^ -427...)))) . (3.23) 

Observe that this method for estabhshing (3.14) avoids the necessity of proving 
hniK^oo Fk = in (3.12) because the right side of (3.12) vanishes exphcitly in this hmit. 



Special Case Rn=n 

An interesting special and simple case of (3.14) arises if we set Rn = n.^ Note that 
Rn = n corresponds to equally-spaced concentric spheres. For this case 

= ^ - m^) • 

where ( is the Riemann zeta function. For this case we have established that for D > 2, 
Fk vanishes as K — > 00. Specifically, we find that for Rn — n, Fk vanishes like K'^~^ as 

The structure of (3.24) is reminiscent of an interesting probabilistic result in number 
theory. It is not difficult to show that the probability of a randomly chosen integer being 
evenly divisible by (contains as a factor) the A'"th power of some integer is 1 — 1/C,{N) for 
N >2. 

Although the solution to the partial difi'erence equation in (3.3) is complicated in 
general for arbitrary D, when R^ = n we have been able to obtain exact solutions for the 
special cases of integer dimension D = 0, 1 and 2. The simplest case is D = 1. For this 
case 

(n + 2j-l)! 



t!t!22* 

Using the generating function formula 



Cn,n+2i-l- ^^^^._^^,^.,2„+2,-2 > 2) 

r _ (2*)! 



(3.25) 



t=i 

which corresponds to (2.10), we find that 



Using (3.27) we verify directly that 



t=l 



\x=l 



1 dx ^ T(t+l/2)x'^' 



Jo \VT^ J 



1 
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For the case = 2, we have 



{D = 2) 



(2n- l)(n + 2j - 1)! 
j!(2n + 2j - l)mi 



Thus, 



C 



1 



l,2t 



l + 2t 



(3.29) 



(3.30) 



Using (3.26) we can obtain from (3.30) the probabihties P2t and we find that they are not 
simple numbers: 

1 

3 ' 
4 

45 ' 
44 



P2 
P4 

Pe 
Ps 
Pio 

Pl2 

Nevertheless, we easily calculate 112: 

X 4x^ 



945 ' 
428 

14175 ' 
10196 

467775 ' 
10719068 

638512875 



n. = - + 

= 1 - 



+ 



.?> 



45 945 



+ 



-^"•"S"*" 5 7 9 



x=l 



1 - 



= 1 . 



y=y/x=l 



(3.31) 



(3.32) 



The case D = is slightly more complicated. For this case there is no simple formula 
for Cn^t- Rather, we must consider special cases: For n > 1 we have 

Cn,n-1 = Bn,0 ; 

Cn,n+1 = Bn,l + g-Bn,0 , 

2 26 (3.33) 

<^n,n+5 — -Dn,3 + ''^ 45 '"^ 945 '° ' 
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and so on, where 

_ (n + 2i-2)!(2n-l 

-DnJ — 



j!(2n + 2j -l)!!2i 
The numbers Cn,t corresponding to n = 1 are 

Ci,o = 1 , 
2 



Ci^2 
(^1,4 



3 ' 
26 

45 ' 

502 (3.34) 

945 ' 
7102 



14175 
44834 



93555 

In contrast, the formula for P2t is extremely simple: 



(2t+l)(2t-l) 

Hence, 



= 77^— -tL^-tt • (3.35) 



Ho = V 7 77 7 = 1 • (3.36) 

^ (2t + l)(2t-l) ^ ^ 



Reciprocity Property of (3.3) 

Finally, we mention an interesting reciprocity property that is exhibited by the prob- 
abilities Ci^2t and P2t- This property is associated with reflecting the dimension D about 
the point D — 1: 



2t 



D 



l,2t-2 



Ci 



,2t\ 



2-D 



(t>l) . 



(3.37) 



An immediate consequence of this probability can be obtained by summing (3.37) from 
t = 1 to oo: 



D 



2-D 



t=l 

= [C'l.o — C'l.c 
= 1-^1, 

Hence Hd = 1 if and only if lim^^oo (^1,2* = in (2 — D) dimensions. 
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(3.38) 



2-D 



We can prove the reciprocity relation in (3.37) in completely general terms. Let Cn,t 
be the solution of (3.3-4), and let C^^t be the solution of (3.3-4) with Pin(n) and Pont{n) 
interchanged for each n > 2 but not n = 1. Then, X^t^o ^i,2t^* is given by (3.19) with 
defined in (3.18), and ^^0^1.2*^* is given by (3.19) with Qn replaced by 



Qn = Pont{n + l)PUn) (n>2) , 

Ql = -Pout (2) . 

With these definitions we obtain according to an identity given by Wall:-'^-'^ 
^ =1/(1 - Q^x/{1 - Q2x/{1 - Q^x/{1 - Q^xj . . .)))) 



(3.39) 



1-a; " (3.40) 

X 1/(1 - Q^xjil - Q2x/{1 - Q^x/{1 - Q^xj ...)))) . 

From (3.19) we then get 

^=(|E^i,2.x*) • (3-41) 

Using the generating function in (3.26) we obtain 

00 00 

Y,P2tx' -\ = {x-\)Y,Cxatx' ; (3.42) 

the coefficient of in this equation is 

P2t = Cx,it-i-Cx,2t (i>l) • (3.43) 

For the particular case in (3.2), interchanging P\yJ^n) and Po^tiji) for n > 2 corresponds to 
replacing D by (2 — D) and relation (3.37) follows immediately. 
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APPENDIX 



In Sec. II we found an expression in (2.12) for the probability IId of a random walker 
eventually returning to the origin on a hypercubic lattice in arbitrary integer dimension D. 
For real values of D > 2, (2.12) describes an acceptable analytic continuation off the 
integers for such a probability. But at D — 2, (2.12) approaches the correct result 112 = 1 
from above with a nonzero slope of — 7r/2. Thus, the physical probability must have a 
discontinuous derivative as a function oi D aX D = 2. Here we calculate an asymptotic 
series for the function that analytically continues (2.12) smoothly to its nonphysical values 
for D <2. The results indicate that this function has a singularity at D = 1. Furthermore, 
this asymptotic series provides excellent numerical results for all D > 1 which can be 
systematically improved at least for integer values of D. 

We start by rewriting the integral representation of Co,2t in (2-7) (recall that Co,t = 
when t is odd): 

/i " 
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2tln 
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Now we consider Co,2t in the limit of a large number of time steps, t — > oo. In this limit 
we can evaluate the integral in (A.l) by Laplace's method^^. The exponent in the integral 
has a maximum either when all D components of k are small or when all D components 
of k are close to tt. Both of these regions give identical contributions to the integral, and 
we can write 
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We scale each of the integration variables in (A. 2) by \/t/D and collect terms order by 
order in powers of 1/t. Note that the choice of e in the end allows one to extend the 
integration region out to infinity to within exponentially small corrections: 
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The remaining integrals separate into products of ordinary integrals to any order in powers 
of 1/t and can be solved with simple combinatorical considerations. We finally obtain 



/ n \ -0/2 oo 



{AA) 



where the Wn{D) are polynomials in D having rational coefficients: 

woiD) = 1 , 

wiiD) = -^D , 

^ ^ 96 128 192 
''^ ^ 256 1,024 512 

' 3,840 18,432 12,288 294,912 40,960 73,728 
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691 ^ 69,760,249 38,695,325,513 . 

wi2(D) = D H D-^ 

^ ' 268,369,920 2,288,988,979,200 1,245,708,288,000 

_ 482,851,775,889,419 ^4 3,429,357, 841,822,897 ^5 

~ 376,702,186,291,200 ^ 273,965,226,393,600 

1, 001, 400, 307, 492, 157, 137^6 2, 388, 031, 903, 466, 391, 071^7 

19,725,496,300,339,200 ^ 21,917,218,111,488,000 

_ 20, 467, 165, 251, 552, 205, 481^8 7, 793, 853, 343, 774, 704, 271^9 

147, 283, 705, 709, 199, 360 ^ 70,135,097,956,761,600 

_ 94, 989, 731, 986, 215, 946, 793 99, 366, 320, 397, 537, 803, 389 ^ 

~ 1,683,242,350,962,278,400 ^ 5,610,807, 836, 540,928,000 

_ 2, 787, 056, 465, 380, 319, 491, 493^^2 10, 961, 446, 588, 181, 276, 321^^3 

~ 888, 751, 961, 308, 082, 995, 200 ^ 45, 848, 315, 464, 305, 868, 800 

24,950,506,464,410,431 292,872,847,633 

^ 103, 687, 728, 819, 276, 349, 440 ^ 1,402,701,959,135,232,000 

318,805,229 59 

1, 683, 242, 350, 962, 278, 400 445, 302, 209, 249, 280 

+ ^ D^"" , 

36,069,478,949,191,680 

(A.5) 

and so on. The polynomials Wn{D) have a few remarkable properties. For instance, the 
degree of the nth polynomial is [3n/2]. The coefficient of D in the nth polynomial, n > 1, 
is Bn2~'^~^ /n where the are the Bernoulli numbers. Also, the coefficient of D^n/s -^-^ 
the 2nth polynomial is 192~'^/n! and the coefficient of Di+3"^/2 ^j^g 2n + 1st polynomial 
is %^192-'^. 

8n! 

To obtain an approximation for Hd in (2.12) we sum (A. 4) over t. Using the definition 
of a generalized Riemann ^-function, 

oo 

as,q) = J2n-', (A6) 

n=q 

we define for N >0: 

HE{x) ^ 2(-] Yl «'-(^) + N+l)x- . (A7) 

^ n=0 

Then we expect that 

~ 1 n (N^oo) . (A.8) 

We obtain excellent numerical results from (A.8) by the following procedure based on 
"optimal truncation" For given D and fixed N truncate the series for H^{x) in (A. 7) 
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right before the term with the smaUest coefficient of x"' and evaluate the truncated series 
ai X = 1. Inserting this value for H^{1) in (A. 8) provides a very accurate numerical 
approximation for IId. The relative size of the first omitted coefficient in the series in 
(A. 7) usually gives a good estimate of the relative error for the value of n^. Fig. 3 
contains the plot of Ho in (A. 8) for A?" = 3 as a function of D for 1 < D < 2. The choice 
of AT = 3 ensures that the error in this plot is well below the width of the line. Note the 
smooth continuation at D = 2 with the plot ofUn in (2.12). 

The function in (A. 8) appears to have a singularity at D = 1: For increasing A^, the 
sum of the two terms X^^=o ^o,2t and (1) in (A. 8) appears to approach zero. This 
cancellation is nontrivial because at D = 1 the former is a rational number while the latter 
is a sum over ^-functions of half-integer arguments. For N = 1, the cancellation is already 
accurate to 0.002%. For N = 15, for example, we find that 

15 

J2 Co,2t = 300, 540, 195/67, 108, 864 
t=o 

and when this rational number is rewritten in decimal form the first 16 digits are canceled 
by Hl^{l). A numerical study indicates that the singularity is a simple pole with a residue 
close to 2. 

For integer values of D, this approximation can be improved systematically by in- 
creasing A^. As an example, we have calculated the value of Us in (2.14) to 0.000003% 
accuracy for N = A. In Table 2 we list the values thus obtained for lis for A/" = 1, 2, 3 and 
4 and the order at which the truncation occurs. Note that even N = 1 gives an accuracy 
of 0.01% in this case. In Tabic 3 we list the values of Ud for the integer dimensions from 
D = 3 to 10 that were obtained up to the indicated accuracy with this method. 

For noninteger values of D, increasing A^ generally does not improve the numerical 
accuracy for in (A. 8). For a given noninteger value of D, N eventually becomes so 
large that Co,2N ceases to be a probability, as we pointed out in Sec. II. Apparently, 
the asymptotic series in (A. 7) cannot correct for this effect and the values for Hd finally 
diverge. 
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such a square is 1/4 x 1/2 = 1/8. This difference reflects the fact that the hypercubic 
lattice is not spherically symmetric. In the spherical geometry pictured in Fig. 5 the 
random walker is equally likely to be at any point in the volume between the spheres. 

Spherical surface areas in this model of random walks determine the relative prob- 
ability of going outward or inward. However, this is not the only way to construct 
a model of a random walk. For an alternative random walk we could generalize the 
notion of Buffon's needle: Suppose we are in region n. We randomly drop a needle 
whose length equals — Rn-i, the radial extent of the region, in such a way that 
one end of the line segment lies in region n. Then, the random walker enters the 
region in which the other end of the line segment lies. For very large concentric circles 
(where the radius is large compared with the spacing) the circumference is locally 
flat and the probability of moving outward, inward, or remaining in the same region 
approaches the usual Buffon's needle result that is discussed in elementary calculus 
textbooks such as G. Klambauer, Aspects of Calculus (Springer, New York, 1986), p. 
365. 
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nonlinear transforms. See C. M. Bender and K. A. Milton, J. Math. Phys. (to be 
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Here, we regard this case as special because in some sense it represents a L>-dimensional 
random walk in Rat space. We can think of curved space as having the property that 
the surface area of the consecutive nested spheres does not grow like n^~^. 
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Entries in a table of values for the probabilities Cn,t form an infinite triangle on an 
(n, t) grid. For a given value of n > 1, n — 1 is the smallest value of t for which Cn.t 
is nonzero. (This is because the random walker cannot reach the nth region until 
time n — 1.) Thus, the lowest value of t contributing to the sum in (3.7) is n — 1. If 
we examine this sum numerically we find that the terms in the sum at first increase 
with t and reach a maximum and then die off to as t continues to increase. Also, 
for fixed t as we vary n, we again observe a local maximum. Thus, in the table of 
probabilities there is a crest; this crest lies on a curve that has the shape of a parabola: 
t{D — 1) = n^. This parabola is just the curve that corresponds to similarity solutions 
of the continuous version (the parabolic partial differential equation) of (3.3). 

Relation (3.40) can be found in H. S. Wall, Analytic Theory of Continued Fractions 
(Van Nostrand, New York, 1948), p. 281, where we identify gn-i = Pinin), n > 2. 

Elementary presentation of Laplace's method may be found in C. M. Bender and 
S. A. Orszag, Advanced Mathematical Methods for Scientists and Engineers (McGraw- 
Hill, New York, 1978), Chap. 6. 

See Ref. 12, Chap. 3, for a discussion of optimal truncation of asymptotic series. 
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FIGURE CAPTIONS 



Figure 1. The first three nonzero probabihties Co,t of returning to the origin after t = 
2, 4, and 6 time steps, plotted as functions of D. Observe that Co,t is not an acceptable 
probability until D is sufficiently large. 

Figure 2. Same as in Fig. 1 except that we have plotted Co, is and Co, 20- For small D 
these functions become so large that it is necessary to scale the vertical axis using the 
transformation y — > [2y/(l + \y\)]. That is, C = ±00 corresponds to ±2 on the vertical 
axis while the interval between and 1, where C is a probability, is mapped onto itself. 

Figure 3. A plot of Hd in (2.12) representing the probability that a random walker on 
a D-dimensional hypercubic lattice will eventually return to the origin. The curve passes 
through 1 at = 2 with a slope of —n/2 and for large D it falls off like 1/{2D). For 
< -D < 2, Ud = 1. However, if we analytically continue the function in (2.12) to values 
of D less than 2 we obtain a function that grows as D decreases and eventually blows up 
at D = 1 (see Appendix). 

Figure 4. Dual hypercubic lattice for random walks in two-dimensional space. A random 
walker begins in the central square labeled 1. On her first step she enters one of the four 
squares labeled 2. On her next step the walker either returns to square 1 or enters one 
of the eight squares labeled 3. In general, at each step the walker goes to an adjacent 
region, either inside or outside her current region. This way of visualizing a random walk 
is completely equivalent to the conventional way in which the walker goes from vertex to 
vertex on a hypercubic lattice. 

Figure 5. The spherical geometry for the arbitrary-dimensional random walks discussed 
in this paper. The random walker begins in the central region and on her first step enters 
region 2. On her second step she either returns to region 1 or moves outward to region 3 
with a relative probability equal to the surface area between her current region and the 
adjacent region to which she goes. 

Figure 6. The first three nonzero probabilities Ci,t of returning to the origin after t = 
2, 4, and 6 time steps, plotted as functions of D, for a random walk in a geometry of 
equally-spaced concentric spheres. Observe that in contrast with Figs. 1 and 2, Ci,t is an 
acceptable probability for all values of D. 

Figure 7. A plot of Hr) in (3.24) representing the probability that a random walker on a 
D-dimensional concentric spherical geometry with equally-spaced spheres will eventually 
return to the origin. Note that the curve is similar to that in Fig. 3; it passes through 1 
at L> = 2 with a slope of -1 and for large D it falls off like 2^"^. For < < 2, = 1. 
If we analytically continue the function in (3.24) to values of D less than 2 we obtain a 
function that grows as D decreases but remains finite until D = —1. AtZ) = lit passes 
through 3 and at = it passes through 13. There are simple poles for all negative 
odd- integer values of D. 
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Table 1 



Table 1. The critical value of the dimension, -Dcriticah for which Co,t can first be considered 
a probability; that is, for all D larger than -Dcriticai, C'o,* lies between and 1. 



critical 



2 0.500 

6 0.722 

10 1.285 

14 1.795 

18 1.990 

22 2.444 

26 3.488 

30 3.892 

34 3.997 

38 4.000 

42 5.647 

46 5.955 

50 6.000 

54 7.243 

58 7.777 
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Table 2 

Table 2. Convergence of the "optimal truncation" procedure to calculate Hd for D = 3. 
The exact result, 0.34053732 . . ., is given in (2.14). 



truncated 

N 113 after order 



1 0.340492587 4 

2 0.340536478 7 

3 0.340537202 9 

4 0.340537315 11 
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Table 3 

Table 3. Values of in (2.12) for some integer values of D. 



D Hd 



3 0.340537329551 

4 0.19320167322 

5 0.135178609 

6 0.104715496 

7 0.085844934 

8 0.07291265 

9 0.06344775 

10 0.05619754 
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